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Abstract 

While the shell model Monte Carlo approach has been successful in the microscopic calculation of nuclear state 
densities, it has been difficult to calculate accurately state densities of odd-even heavy nuclei. This is because the 
projection on an odd number of particles in the shell model Monte Carlo method leads to a sign problem at low 
temperatures, making it impractical to extract the ground- state energy in direct Monte Carlo calculations. We show 
that the ground-state energy can be extracted to a good precision by using level counting data at low excitation 
energies and the neutron resonance data at the neutron threshold energy. This allows us to extend recent applications 
of the shell-model Monte Carlo method in even-even rare-earth nuclei to the odd-even isotopic chains of 9 ~ 55 Sm 
and I43 ~ I49 Nd. We calculate the state densities of the odd-even samarium and neodymium isotopes and find close 
agreement with the state densities extracted from experimental data. 

Keywords: nuclear shell model, quantum Monte Carlo, level density, back-shifted Bethe formula 
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1. Introduction 

Reliable microscopic calculation of nuclear state den- 
sities of heavy nuclei is a challenging task because it 
often requires the inclusion of correlations beyond the 
mean-field approximation. Correlation effects can be 
taken into account in the context of the configuration- 
interaction (CI) shell model approach, but the size of the 
required model space in heavy nuclei is prohibitively 
large for direct diagonalization of the CI shell model 
Hamiltonian. This limitation can be overcome in part by 
usingthe shell model Monte Carlo (SMMC) method fl 
SSSSl. The SMMC method enables the calculations 
of statistical nuclear properties, and inparticular of level 
densities in very large model spaces lid. MM foil id]. 

Fermionic Monte Carlo methods are often hampered 
by the so-called sign problem, which leads to large and 
uncontrollable fluctuations of observables at low tem- 
peratures yj, |2l [ill UM- Statistical and collective prop- 
erties of nuclei can be reliably calculated [ 13, 6] by em- 
ploying a class of interactions that have a good Monte 
Carlo sign in the grand-canonical formulation lul Il3ll. 
The projection on an even number of particles keeps 
the good sign of the interaction, allowing accurate cal- 



culations for even-even nuclei. However, the projec- 
tion on an odd number of particles leads to a new sign 
problem, making it difficult to calculate thermal observ- 
ables at low temperatures for odd-even and odd-odd nu- 
clei. In particular, the ground-state energy of the odd- 
particle system cannot be extracted from direct SMMC 
calculations. A method to extract the ground-state en- 
ergy of the odd-particle system from the imaginary-time 
Green's functions of the even-particle system was re- 
cently introduced in Ref. [ 14] and applied successfully 
to medium-mass nuclei. However, the application of 
this method to heavy nuclei is computationally inten- 
sive and requires additional development. In this work, 
we describe a practical method that enables us to de- 
termine the ground-state energy using a complete set 
of experimentally known low-lying levels and the neu- 
tron resonance data at the neutron threshold energy. We 
then extend the recent SMMC state density calculations 
of even-even nuclei in the rare-earth region 11151 1 1611 
to include odd-even isotope chains of 149 ~ I55 Sm and 
143 ~ I49 Nd. We find close agreement with the state den- 
sities that are extracted from experimental data. 
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2. Choice of Model Space and Interaction 



3. Ground-state Energy 



In rare-earth nuclei we use the model space of 
Refs. I115[ llql spanned by the single-particle orbitals 
0g 7 /2, ld 5 /2, ldi/2, 2*1/2, O/111/2 and l/ 7/2 for protons, 
and the orbitals 0h n/2 , 0h 9/2 , I/7/2, I/5/2, 2p 3/2 , 2p l/2 , 
O/13/2, and lg9/2 for neutrons. We have chosen these 
orbitals by the requirement that their occupation proba- 
bilities in well-deformed nuclei be between 0.9 and 0. 1 . 
The effect of other orbitals is accounted for by the renor- 
malization of the interaction. The bare single-particle 
energies in the shell-model Hamiltonian are obtained so 
as to reproduce the single-particle energies in a spheri- 
cal Woods-Saxon plus spin-orbit potential in the spher- 
ical Hartree-Fock approximation. 

The effective interaction consists of monopole pairing 
and multipole-multipole terms 111511 
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where the pair creation operator P v and the multipole 
operator Oa- v are given by 
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with CI; 



(2b) 
*jm - {—y +m dj- m - In Eq. ([1]), : : denotes normal 
ordering and Vws represents the Woods-Saxon poten- 
tial. The pairing strengths is expressed as g v - yg v , 
where y is a renormalization factor and g p - 10.9/Z, 
g„ = 10.9 /N are parametrized to reproduce the exper- 
imental odd-even mass differences for nearby spher- 
ical nuclei in the number-projected BCS approxima- 
tion 111511 . The quadrupole, octupole and hexadecupole 
interaction terms have strengths given by xa - &AX f° r 
A = 2,3,4 respectively. The paprameter x is deter- 



mined self-consistently IU3I1 and k* are renormalization 
factors accounting for core polarization effects. We use 
the parametrization of Ref. [16] for y and k 2 
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while the other parameters are fixed at £3 = 1 and k\ 
1. 



Because of the sign problem introduced by the pro- 
jection on an odd number of particles, the thermal en- 
ergy E(J3) of the odd-even samarium and neodymium 
isotopes can in practice be calculated only up to /3 = 
1/r ~ 4-5 MeV _1 . For comparison, SMMC cal- 
culations in the neighboring even-even samarium and 
neodymium nuclei, for which there is no sign problem, 
were carried out up to (3 — 20 MeV -1 B16I1 . System- 
atic errors introduced by the discretization of /3 yfl, are 
corrected by calculating E(J3) for the two time slices of 
A/3 = 1/32 MeV" 1 and A/3 = 1/64 MeV" 1 and then per- 
forming a linear extrapolation to A/3 - 0. For /? ~ 3 
MeV~', the dependence of E(J3) on A/3 is weaker and 
an average value was taken instead. The calculations 
for jS > 2 MeV~' were carried out using a stabilization 
method of the one-body canonical propagator [15] for a 
given configuration of the auxiliary fields. 

Since our calculations of the thermal energy E(J3) are 
limited to/3 ~ 4 — 5 MeV -1 , we cannot obtained a re- 
liable estimate of the ground-state energy Eq in direct 
SMMC calculations. The Green's function method of 
Ref. 114(1 was used successfully in medium-mass nuclei 
to circumvent the odd particle-number sign problem and 
extract accurate ground-state energies. However, this 
method becomes computationally intensive in heavy nu- 
clei and requires additional development. Here we ex- 
tract Eq by performing a one-parameter fit of the SMMC 
thermal excitation energy E X {T) = E(T) - Eq to the ex- 
perimental thermal energy. The latter is calculated using 
the thermodynamical relation E X (J3) = —d \n Z(J3)/d/3, 
where Z is the experimental partition function in which 
the energy is measured relative to the ground state (see 
below). 

We demonstrate the method for I47 Nd. Fig. [T] shows 
the thermal excitation energy (left panel) and the par- 
tition function (right panel) as a function of temper- 
ature. The SMMC results (open circles) are shown 
down to T - 0.25 MeV, below which the statistical er- 
rors become too large due to the odd particle-number 
sign problem. We calculated the experimental parti- 
tion function (right panel, dashed line) from Z(T) - 
2,(2-/; + l)e~ Ex -^ T where E x j denote the excitation en- 
ergies of the experimentally known energy levels. The 
incompleteness of the level counting data above a cer- 
tain excitation energy leads to a saturation of the ex- 
perimental partition function and the corresponding av- 
erage experimental thermal energy (left panel, dashed 
line) above T ~ 0.15 MeV. In order to determine the 
ground-state energy Eq from a fit of the SMMC ther- 
mal energy to the experimental data, it is necessary to 
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Figure 1 : Thermal excitation energy E x (left panel) and partition func- 
tion Z (right panel) versus temperature T for l47 Nd. The SMMC re- 
sults (open circles) are compared with the results deduced directly 
from experimentally known levels (dashed lines), and from Eq. (4} 
(solid lines). The ground-state energy Eq is obtained by a one- 
parameter fit of the SMMC thermal energy to the solid line on the 
left panel. 



have a realistic estimate of the experimental thermal en- 
ergy above T ~ 0.25 MeV where SMMC results are 
available. This can be accomplished by calculating an 
empirical partition function 

Z(T) = YpJi + V)e~ E " IT + \ dE xPBB¥ (E x )e- E/T , 

(4) 
where the discrete sum is carried out over a suitably cho- 
sen complete set of N experimental levels up to an ex- 
citation energy En, and the contribution of higher-lying 
levels is included effectively through the integral over 
an empirical level density with a Boltzmann weight. 
The empirical state density is parametrized by the back- 
shifted Bethe formula (BBF) Q3] 
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where a is the single-particle level density parameter 
and A is the backshift parameter. For these parameters 
we use the values in Ref. 11811 determined from level 
counting data (at low excitation energies) and the s- 
wave neutron resonance data (at the neutron resonance 
threshold) for the given nucleus of interest. The solid 
lines in Fig. Q] are calculated using Eq. (|4]i for the em- 
pirical partition function. The number of the low-lying 
states N (determining En) are chosen such that the solid 
and dashed curves for the experimental partition func- 
tion merge smoothly at sufficiently low temperatures. 
The values of a and A and the suitably chosen values 
of N are tabulated in Table Q] for the nuclei of interest. 
The ground-state energy Eq is determined by fitting the 
SMMC thermal excitation energy (open circles) to the 
solid curve in the left panel of Fig. Q] TableQ~]also shows 
the fitted values of Eq for the corresponding neodymium 
and samarium isotopes. 
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Figure 2: State densities of odd-mass 143 ~ 149 Nd isotopes versus exci- 
tation energy E x . The SMMC densities (open circles) are compared 
with experimental results. The histograms at low excitation energies 
are from level counting data, the triangle is the neutron resonance data 
and the dashed line is the BBF state density [Eq. J3J] that is extracted 
from the experimental data. 



We emphasize that only a single parameter Eq is ad- 
justed to fit the empirical partition function. Neverthe- 
less, we see that both the SMMC partition function and 
thermal energy fit very well the corresponding experi- 
mental quantities over a range of temperatures. This is 
an evidence that our SMMC results are consistent with 
the experimental data. 

4. State Densities 

The average state density is determined in the saddle- 
point approximation to the integral that expresses the 
state density as an inverse Laplace transform of the par- 
tition function. This average density is given by 
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where S(E) is the entropy and C is the heat capacity 
in the canonical ensemble. In SMMC we first calculate 
the thermal energy as an observable E(J3) - (H) and the 
integrate the thermodynamic relation -dlnZ/d/3 - E(J3) 
to find the partition function Z(J3). The entropy and the 
heat capacity are then calculated from 



S (E) = InZ + pE ; C = ~p 2 dE/dp 



(7) 



Nucleus 


a (MeV- 1 ) 


A (MeV) 


N 


E N (MeV) 


E (MeV) 


143 Nd 


15.951 


0.759 


3 


1.228 


-191.609 ±0.020 


145 Nd 


16.792 


0.135 


4 


0.658 


-210.010 ±0.028 


147 Nd 


18.276 


0.058 


4 


0.190 


-227.870 ± 0.029 


149 Nd 


18.552 


-0.644 


3 


0.138 


-245.440 ±0.021 


149 Sm 


19.086 


-0.196 


8 


0.558 


-244.894 ± 0.026 


15 'Sm 


18.999 


-0.771 


8 


0.168 


-264.968 ± 0.020 


153 Sm 


17.848 


-1.053 


7 


0.098 


-284.796 ± 0.020 


155 Sm 


17.067 


-0.834 


7 


0.221 


-304.245 ± 0.025 



Table 1 : The values of a and A in the BBF [see Eq. (3)] as determined from the level counting data at low excitation energies and the neutron 
resonance data 1 18], the number N of a complete set of experimentally known levels and its corresponding energy En used in Eq. {4), and the 



extracted ground-state energies Eq for the odd-mass ; 



y Nd and ' 



Sm isotopes. 



and substituted into Eq. (O to yield the state den- 
sity. The SMMC state densities for the odd-mass 
neodymium and samarium isotopes are shown by open 
circles in Fig. [2] and Fig. [3] respectively. We com- 
pare these SMMC densities with the level counting data 
(histograms) and with the neutron resonance data (tri- 
angles). We also show the empirical BBF level densi- 
ties (dashed lines). For the neodymium nuclei (Fig. [2j, 
we find excellent agreement among the SMMC results 
(open circles), the experimental state densities (both at 
the level counting and at the neutron resonance ener- 
gies) and the BBF state densities (dashed lines). In 
the case of the samarium nuclei (Fig. [3}, a similarly 
good agreement is observed for 149 Sm and 151 Sm. We 
also find overall good agreement for 153 Sm and 155 Sm, 
although we observe some discrepancies between the 
SMMC and the experimental state densities at the neu- 
tron resonance energy. Since the neutron resonance data 
point is used to determine the parameters of the BBF 
state densities, similar discrepancies are observed be- 
tween the SMMC and the BBF state densities. 
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Figure 3: State densities of odd-mass ' 155 Sm isotopes versus exci- 
tation energy E x . Symbols and lines as in Fig. [2] 



5. Conclusion 

We have extended the SMMC calculations in the rare- 
earth region to include the odd-mass isotopic chains of 
i49-i55o; m an( j 143-149^ -pj^ s jg n p ro blem introduced 

by the projection on an odd number of particles makes 
it difficult to calculate directly the ground-state energy. 
We circumvent this problem in practice by carrying out 
a one-parameter fit of the SMMC thermal excitation en- 
ergy to the experimental average thermal excitation en- 
ergy as determined from the neutron resonance data and 
level counting data at low excitation energies. We then 
calculate the state densities of the even-odd samarium 
and neodymium isotopes and find them to be in good 



overall agreement with state densities extracted from 
available experimental data. This method requires suffi- 
cient level counting data at low excitation energies and 
the neutron resonance data point at the neutron reso- 
nance threshold. The good overall agreement between 
the SMMC and experimental densities is an evidence 
that the theory is consistent with the data. 
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